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1. Intrinsic heavy flavours in the proton 

The Large Hadron Collider (LHC) at CERN can be a useful laboratory for investigation of the 
unique structure of the proton, in particular for the study of the parton distribution functions (PDFs) 
with high accuracy. It is well known that the precise knowledge of these PDFs is very important 
for verification of the Standard Model and search for New Physics. 

By definition, the PDF fa{x,}x) is a function of the proton momentum fraction x carried by 
parton a (quark q or gluon g) at the QCD momentum transfer scale /i. For small values of ji, 
corresponding to the long distance scales less than I//X0, the PDF cannot be calculated from the 
first principles of QCD (although some progress in this direction has been recently achieved within 
the lattice methods [[T]]). The PDF /q(x,/x) at > /Iq can be calculated by solving the perturbative 
QCD evolution equations (DGLAP) The unknown (input for the evolution) functions /a(x,/Xo) 
can usually be found empirically from some "QCD global analysis" ^] of a large variety of data, 
typically at /i > /lo- 
in general, almost all pp processes that took place at the LHC energies, including the Higgs 
boson production, are sensitive to the charm /c(x,/i) or bottom /fo(x,/i) PDFs. Nevertheless, within 
the global analysis the charm content of the proton at /i ~ jic and the bottom one at ~ /x^ are 
both assumed to be negligible. Here pLc and pL}, are typical energy scales relevant to the c- and 
Z?-quark QCD excitation in the proton. These heavy quark components arise in the proton only 
perturbatively with increasing 2^-scale through the gluon splitting in the DGLAP evolution 
Direct measurement of the open charm and open bottom production in the deep inelastic processes 
(DIS) confirms the perturbative origin of heavy quark flavours [||]. However, the description of 
these experimental data is not sensitive to the heavy quark distributions at relatively large x (x > 
0.1). 

As was assumed by Brodsky with coauthors in 0], there are extrinsic and intrinsic con- 
tributions to the quark-gluon structure of the proton. Extrinsic (or ordinary) quarks and gluons 
are generated on a short time scale associated with a large-transverse-momentum processes. Their 
distribution functions satisfy the standard QCD evolution equations. Intrinsic quarks and gluons 
exist over a time scale which is independent of any probe momentum transfer. They can be asso- 
ciated with bound-state (zero-momentum transfer regime) hadron dynamics and are believed to be 
of nonperturarbative origin. Figure |T] gives a schematic view of a nucleon, which consists of three 
valence quarks qy, quark-antiquark qq and gluon sea, and, for example, pairs of the intrinsic charm 
(lin'lin) ^'^'^ intrinsic bottom quarks (qS'^qj't,)- 

It was shown in [^] that the existence of intrinsic heavy quark pairs cc and bb within the proton 
state could be due to the virtue of gluon-exchange and vacuum-polarization graphs. On this basis, 
within the MIT bag model [^, the probability to find the five-quark component \uudcc) bound 
within the nucleon bag was estimated to be about 1-2%. 

Initially in ||^ ^ S. Brodsky with coauthors have proposed existence of the 5-quark state 
\uudcc) in the proton (Fig. [l]). Later some other models were developed. One of them consid- 
ered a quasi-two-body state D^{uc)A^{udc) in the proton [^. In [@l-[|ll|] the probability to find 
the intrinsic charm (IC) in the proton (the weight of the relevant Fock state in the proton) was as- 
sumed to be 1-3.5%. The probability of the intrinsic bottom (IB) in the proton is suppressed by 



the factor m^/m^ ~ 0.1 [12], where nic and nib are the masses of the charmed and bottom quarks. 
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Figure 1: Schematic presentation of a nucleon consisting of three valence quarks qy, quark-antiquark qq 
and gluon sea, and pairs of the intrinsic charm (qf.,qf„) and intrinsic bottom quarks (q|lqf!,)- 



Nevertheless, it was shown that the IC could result in a sizable contribution to the forward charmed 



meson production [13]. Furthermore the IC "signal" can constitute almost 100% of the inclusive 



spectrum of D-mesons produced at high pseudorapidities Tj and large transverse momenta pr in 



pp collisions at LHC energies [14] 



If the distributions of the intrinsic charm or bottom in the proton are hard enough and are 
similar in the shape to the valence quark distributions (have the valence-like form), then the pro- 
duction of the charmed (bottom) mesons or charmed (bottom) baryons in the fragmentation region 
should be similar to the production of pions or nucleons. However, the yield of this production 
depends on the probability to find the intrinsic charm or bottom in the proton, but this yield looks 
too small. The PDF which included the IC contribution in the proton have already been used in the 
perturbative QCD calculations in [^-[|n]]. 

The probability distribution for the 5-quark state (uudcc) in the light-cone description of the 
proton was first calculated in Assuming that the light quark {u,d) masses and the proton 



mass are smaller than the c-quark mass one can get the following form for this probability [IC] at 
= (nic = 1.69 is the c-quark mass): 

dP 

— =f,{x)=Mx)=J^x^{{l-x){l + \0x + x^) + 6x{l+x)ln{x)} , (1.1) 

where the normalization constant ^ determines some probability Wj^ to find the Fock state \uudcc) 
in the proton. The solid line in Fig. || shows the intrinsic charm PDF xfc{x) as a function of x at 
= ml when the probability w^q — 3-5%. The dashed curve in Fig. ^ is the density distribution 
of the (ordinary) sea charm in the proton. One can see from Fig. that the IC distribution (with 



w^Q = 3.5% IpTl]) given by Eq. ( [Til ) has rather visible enhancement at x ~ 0.2 — 0.3 and it is much 
larger (a few orders of magnitude) than the sea (ordinary) charm density distribution in the proton. 
The valence quark density distributions {xfuj), the distributions of the sea charm (Csea) and the 
intrinsic charm {xf\c{x)) at = 100 GeV/c^ are presented in Fig. ^ One can see from Fig. ^. that 
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Figure 2: The distributions of charmed quarks in the proton; the dashed line is the sea charmed quarks c{x), 
the solid curve is the sum of the intrinsic charm ans the sea one c(x) + Cin(x) 




Figure 3: The distributions of the valence {u,d) quarks and sea quarks in the proton as a function of the 
Bjorken variable x. 
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the sea charm and the intrinsic charm distributions are much "smaller" than the valence ones in the 
whole region of x. However, in the hard pp collisions the gluons and sea quarks make the main 
contribution to the inclusive hadron spectra. Therefore, the inclusion of the intrinsic heavy quark 
components in the proton makes sense. 

Due to the nonperturbative intrinsic heavy quark components one can expect some excess 
of the heavy quark PDFs over the ordinary sea quark PDFs at x > 0.1. The "signal" of these 
components can be visible in the observables of the heavy flavour production in semi-inclusive 
ep DIS and inclusive pp collisions at high energies. For example, it was recently shown that 
rather good description of the HERMES data on the xfs{x,Q^) +xfs{x,Q^) at x > 0.1 and = 
2.5 GeV/c^ [ p3| ] could be achieved due to existence of intrinsic strangeness in the proton. Similarly, 
possible existence of the intrinsic charm in the proton can lead to some enhancement in the inclusive 
spectra of the open charm hadrons, in particular D-mesons, produced at the LHC in pp-colhsions 



at high pseudorapidities Tj and large transverse momenta pr [14]. 



2. Intrinsic heavy quarks at the LHC 
• Intrinsic cliarm 

It is known that in the open charm/beauty /j/j-production at large momentum transfer the hard 
QCD interactions of two sea quarks, two gluons and a gluon with a sea quark play the main role. 



According to the model of hard scattering []21|]-[[29|] the relativistic invariant inclusive spectrum of 
the hard process p + p ^ h+X can be related to the elastic parton-parton subprocess i + j ^ i' + f, 
where i,j are the partons (quarks and gluons). This spectrum can be presented in the following 



general form [M]-[27] (see also [BO, 31]): 



= I / ^'^'T / d%r f dx, f dxjMxMfj{xj,kjr)^^^^i^. (2.1) 
dp TiJ J A'"'" A'™" dt KZh 

Here kjj and j are the four-momenta of the partons / or j before and after the elastic parton-parton 
scattering, respectively; kij^kjj are the transverse momenta of the partons / and j; z is the fraction 
of the hadron momentum from the parton momentum; fij is the PDF; and D, y is the fragmentation 
function (FF) of the parton / or j into a hadron h. 

When the transverse momenta of the partons are neglected in comparison with the longitudinal 



momenta, the variables s, t, u and Zh can be presented in the following forms [25]: 



t U X\ X7 

S — XjXjS^ I — Xj — ^ U — Xj — J Zli — — ~t~ — ) y^Zd.Zd J 

Zh Zh Xj 

where 

xi = — = yCot(e/2), X2 = — = ytan(0/2), xt = 2Vtu/s = 2pt/\/s. (2.3) 

Here as usual, s = {pi+ P2)^, t = {pi — p\)^, u = {p2 — p\Y, and p\, pi, p\ are the 4-momenta 
of the colliding protons and the produced hadron h, respectively; is the scattering angle for the 
hadron h in the pp c.m.s. The lower limits of the integration in ( |2. 1[ ) are 

min _ ^rcot(f) ■ _ x,xrtan(f) 

2 — xrtan(^) 2x,- — X7Cot(^) 
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Actually, the parton distribution functions fi{xi^kij) also depend on the four-momentum transfer 
squared that is related to the Mandelstam variables s, f , u for the elastic parton-parton scattering 



2s tu 



s^ + F + u^ 



(2.5) 



Calculating spectra by Eq.(2.1) we used the PDF which includes the IS (and does not include 



it) [11], the FF of the type AKK08 [29] and dOij{s,t)/dt calculated within the LO QCD and 



presented, for example, in []28[]. 

In particular. Fig. ^ shows our estimation of the inclusive yield of single D"-mesons in pp- 
production at = 7 TeV and 10 GeV/c < pj <25 GeV/c as a function of the pseudorapidity rj. 
This estimation is obtained within PYTHIAS, where the PDF CTEQ66 (without the IC, given by 



9000 



10GeV<p^<25 GeV 

CTEQ66C 3,5% IC 
Entries: 781710 

CTEQ66 
Entries: 727285 




6.5 



Figure 4: The D + Dq distributions (with and without intrinsic charm contribution) over the pseudorapidity 

T] for pp {Do+Do) +X at = 7 TeV and 10 < pr < 25 GeV/c [Q. 



the dashed blue curve in Fig. ^ and CTEQ66c (with the IC for Wj^ — 3.5%, given by the sohd red 



curve) were used at = ml =(1.69 GeV)^ Ollp 



One can see from Fig. ^ that there is some enhancement due to the IC contribution at large t] , 
which is due to the above-mentioned enhancement in the IC PDF at x > 0.1, given in Fig. ^. Its 
amount increases with growing pj. For example, due to the IC the spectrum increases by a factor 
of 2 at T] =4.5. A similar effect was predicted in [02|]. 

One can see that the Feynman variable xp of the produced hadron, for example, the D°-meson, 
can be expressed via the variables pj and r\,or Q the hadron scattering angle in the pp c.m.s, 

.,.?^i = ?f-l- = ?^si,AC). (2.6, 

y/s y/s tan0 ^JS 
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At small scattering angles of the produced hadron this formula becomes 

2pT 1 

Xf-^^-. (2.7) 

It is clear that for fixed pr an outgoing hadron must possess a very small 6 or very large Tj in order 
to have large Xf (to follow forward, or backward direction). 

In the fragmentation region (of large xp) the Feynman variable Xf of the produced hadron is 
related to the variable x of the intrinsic charm quark in the proton, and according to the longitudinal 
momentum conservation law, the x/r ~ x (and Xf < x). Therefore, the visible excess of the solid 
(red) histogram over the dashed (blue) one in Fig. § at t] > 3.5 is due to the enhancement of the IC 
distribution (see Fig. |2|) at x > 0.1. 

One expects similar enhancement in the experimental distributions of the open bottom produc- 
tion due to the (hidden) intrinsic bottom (IB) in the proton, which could have the PDF very similar 
to the one given in ([ET|). However, the probability wib to find in the proton the Fock state with 
the IB contribution \uudbb) is about 10 times lower than the IC probability w\c due to the relation 
wib/wic ~ m^/m^, where mj, is the bottom quark mass |l2|]. 

• Intrinsic strangeness 

Let us analyze now how the possible existence of the intrinsic strangeness in the proton can be 
visible in pp collisions. For example, consider the /T^-meson production in the process pp — )• 
+ X. Considering the intrinsic strangeness in the proton [|5l] we calculated the inclusive spec- 



trum EDo /d p of such mesons within the hard scattering model (Eq.(2. 1)), which describes sat 



isfactorily the HERA and HERMES data on the DIS. The FF and the parton cross sections were 



taken from [ |29| , |28| ], respectively, as mentioned above. In Figs. (gj6|) the inclusive -spectra of 
i^^-mesons produced in pp collision at the initial energy Ep =158 GeV are presented at the rapid- 
ity y =1.3 (Fig|5|) and y =1.7 (Fig ^. The solid lines in Figs. (||,^ correspond to our calculation 
ignoring the intrinsic strangeness (IS) in the proton and the dashed curves correspond to the cal- 



culation including the IS with the probability about 2.5%, according to [|15[]. The crosses show the 
ratio of our calculation with the IS and without the IS minus 1 . One can see from Figs. (^,^ right 
axis) that the IS signal can be above 200 % at j = 1.3, pt = Gev/c and slightly smaller, 

than 200 %a.ty= 1.7, pt ~ 2.5 Gev/c. Actually, this is our prediction for the NA61 experiment 
that is now under way at CERN. 

We also calculated the inclusive spectra of K~ and mesons produced in the pp collision 
at LHC energies. Note that the spectra of -mesons, which consist of s- and U- quarks, can give 
us information on the IS at large Xf (xf >0.1) or at the certain region of pj and rj, according to 



Eq.([2.q) and Fig. (|2|), because they are produced mainly from the fragmentation of the strange sea 
(Ssea) and the intrinsic strange (si^) quarks in the proton. One can see from Fig. (^ right axis) that 
the IS signal at ^/s =7 TeV, y =4.5 can be about 400 % at /jy ~ 32-33 GeV/c, the measuring of 
which at LHC can be difficult. Anyway, at pj >15 GeV/c and y =4.5 the IS signal is about 200 % 
and more, as the crosses show in Fig. (fright axis). However, the spectra of ^+-mesons almost do 
not give us such information because ^+ -meson consists of the u- and ^-quarks and its production 
at large xp is due to the fragmentation of the valence w-quark by the pp collision. It illustrates 
Fig. (^, the notations are the same as in Figs. (^^. The similar effect but for the IC signal can be 
visible in the inclusive spectra of D - and D+-mesons produced in pp collision at LHC energies. 
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pp->K"X, Ep=158 GeV, y=1.3 
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Figure 5: The K -meson distributions (with and without intrinsic strangeness contribution) over the trans- 
verse momentum pt for pp K +X at the initial energy E = 158 GeV, the rapidity y =1.3 and pt > 0.8 
GeV/c. 



pp->K" X, Ep=1 58 GeV, y=1 .7 
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Figure 6: The K -meson distributions (with and without intrinsic strangeness contribution) over the trans- 
verse momentum pt for pp K +X at the initial energy Ep = 158 GeV, the rapidity y —1.7 and pt > 0.8 
GeV/c. 



Let us also note that we calculated the inclusive spectra of K- and D-mesons within the LO 
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QCD. The NLO calculation, in principle, can change the inclusive spectra. However, the ratio of 
the spectrum of K (or D ) to the spectrum of (or D+) can not be very sensitive to the NLO 
corrections. 




Figure 7: The K -meson distributions (with and without intrinsic strangeness contribution) over the trans- 
verse momentum p, for pp — > K +X at y/s = 7 TeV, the rapidity y =4.5 and pt > 10 GeV/c. 
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Figure 8: The /r+-meson distributions (with and without intrinsic strangeness contribution) over the trans- 
verse momentum pt for pp — > +X at — 7 TeV, the rapidity y —4.5 and pt > 10 GeV/c. 
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3. Conclusion 



In this paper we have shown that possible existence of the intrinsic heavy flavour quark com- 
ponents in the proton can be seen by the forward production of the open heavy flavour in pp- 
coUisions. 

Our calculations of the charmed meson production in pp collisions within the hard scattering 
model using the PYTHIAS MC generator Q] and the PDF including the intrinsic charm [11] 
showed the foUowing. We found that the contribution of the intrinsic charm in the proton could 
be studied in the production of D-mesons in pp coUisions at the LHC. The IC contribution for the 
single D^'-meson production can be large, about 100 % at high rapidities 3 <y < 4.5 and large 
transverse momenta 10 < < 25 GeV/c. For the double D" production this contribution is not 
larger than 30 % at > 5 GeV/c and 3 <y < 4.5. These IC contributions for the single and 
double D-meson production [|^] were obtained with the probability of the intrinsic charm taken to 
be Wcc =3.5 % [ 1 1 1, and they will decrease by a factor of 3 when Wcs — 1 %. Therefore, this value 
can be verified experimentally at the LHCb. 

The presented predictions could stimulate measurement of the single and double D-meson 
production in pp collisions at the CERN LHCb experiment in the kinematic region mentioned 
above to observe a possible signal for the intrinsic charm. The intrinsic beauty in the proton is 
suppressed by a factor of 10, therefore its signal in the inclusive spectra of B-mesons will probably 
be very weak. 

We also analyzed the inclusive -meson production in pp collision at the initial energy 
Ep =158 GeV and gave some predictions for the NA61 experiment going on at CERN, and made 
some predictions for the inclusive spectra of -mesons produced in pp at the LHC energy ^/s = 
7 TeV. We showed that in the inclusive spectrum of A'^ -mesons as a function of pt at the initial 
energy Ep =15S GeV and some values of their rapidities the signal of the intrinsic strangeness can 
be visible and reach about 200% and more at large momentum transfer we took. The signal of the 
intrinsic strangeness in the inclusive spectrum of A'^ -mesons produced in pp collision at ^/s =7 
TeV can be about 200% and even 400% at pt >15 GeV/c, y =4.5. The probability of the intrinsic 
strangeness to be about 2.5%, as was found from the best description of the HERA and HERMES 
data on the DIS, see [15] and references therein. 
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